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ABSTRACT 

Relativistic and non-relativistic modern nucleon-nucleon potentials are mapped on a rel- 
ativistic operator basis using projection techniques. This allows to compare the various 
potentials at the level of covariant amplitudes were a remarkable agreement is found. In 
nuclear matter large scalar and vector mean fields of several hundred MeV magnitude are 
generated at tree level. This is found to be a model independent feature of the nucleon- 
nucleon interaction. 



1 



I. INTRODUCTION 



A fundamental question in nuclear physics is the role which relativity plays in nuclear 
systems. The ratio of the Fermi momentum over the nucleon mass is about k-p/M ~ 0.25 
and nucleons move thus with maximally about 1/4 of the velocity of light which implies 
only moderate corrections from relativistic kinematics. However, there exists a fundamental 
difference between relativistic and non-relativistic dynamics: a genuine feature of relativistic 
nuclear dynamics is the appearance of large scalar and vector mean fields, each of the 
magnitude of several hundred MeV. The scalar field Eg is attractive and the vector field E M 
is repulsive. In relativistic mean field (RMF) theory, both, sign and size of the fields are 
enforced by the nuclear saturation mechanism 

At nuclear saturation density po — 0.16 fm~ 3 the empirical fields deduced from RMF fits 
to finite nuclei are of the order of S s ~ —350 MeV and S Q ~ +300 MeV 0] (S Q is the time- 
like component of ). The single particle potential in which the nucleons move originates 
from the cancellation of the two contributions U s , p , ~ S D + S s ~ —50 MeV which makes 
it difficult to observe relativistic effects in nuclear systems. There exist, however, several 
features in nuclear structure which can naturally be explained within Dirac phenomenology 
while models based on non-relativistic dynamics have difficulties. Best established is the 
large spin-orbit splitting in finite nuclei. Also the so-called pseudo-spin symmetry, observed 
more than thirty years ago in single-particle levels of spherical nuclei, can naturally be 
understood within RMF theory as a consequence of the coupling to the lower components 
of the Dirac equation ^j]. 

A connection to Quantum-Chromo-Dynamics (QCD) as the fundamental theory of strong 
interactions is established by QCD sum rules 0, 0|. The change of the chiral condensates 
(qq), (q^q) in matter leads to attractive scalar and repulsive vector self-energies which are 
astonishingly close to the empirical values derived from RMF fits to the nuclear chart. Also 
relativistic many-body calculations 0,0,0 yield scalar/vector fields of the same sign and 
magnitude as obtained from RMF theory or, alternatively, from QCD sum rules. Moreover, 
Dirac-Brueckner-Hartree-Fock (DBHF) calculations |2| agree even quantitatively surpris- 
ingly well with the QCD motivated approach of Ref. |9( where chiral fluctuations from 
pion-nucleon dynamics were considered on top of the chiral condensates. 

These facts suggest that preconditions for the existence of large fields in matter or, alter- 
natively, the density dependence of the QCD condensates, must already be inherent in the 
vacuum nucleon-nucleon (NN) interaction. 

Relativistic one-boson-exchange potentials (OBEP), e.g. Bonn (A,B,C), |10J] provide fits 
with fair precision to NN scattering data as well as the non-relativistic meson-theoretical 
Nijmegen soft-core potential Nijm93 [3]. High precision fits (with % 2 /datum ~ 1.01) are 
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FIG. 1: Diagonal matrix elements V(c{, q) in the 1 So partial wave U(q, q) for different high precision 
NN potential models. 



provided by OBE type potentials such as the relativistic CD-Bonn [TJ| or the non-relativistic 
Nijmegen potentials Nijm I and Nijm II [13] where fits are performed separately in each 
partial wave. In contrast to e.g. the spectator model of [14] the 'standard' OBE potentials 
are based on the no-sea approximation which excludes explicit excitations of anti-nucleons 
but absorbs such contributions into the model parameters, in particular into a large u 
coupling. The present investigations are restricted to the standard-type models. 

The long-range part of the interaction is generally mediated by the one-pion-exchange 
(OPE) while the scalar intermediate range attraction is mainly due to correlated two-pion- 
exchange. The short-range part, i.e. the hard core, is dominated by vector meson exchange, 
i.e. the isoscalar uj and the isovector p meson. Widely used in nuclear physics are, however, 
also high precision non-relativistic empirical potentials such as the Argonne potential AVi§ 
(l5l | or the Reid93 potential ^J. A systematic connection to QCD is established by chiral 
effective field theory (EFT). Up to now the two-nucleon system has been considered at 
next-to-next-to-next-to-leading order (N 3 LO) in chiral perturbation theory 0, In such 
approaches the NN potential consists of one-, two- and three-pion exchanges and regularizing 
contact interactions which account for the short-range correlations. The advantage of EFT 
compared to the more phenomenological OBE potentials is the systematic expansion in 
terms of chiral power counting. 

A better understanding of the common features of the various approaches is essential 
in order to arrive at a more model independent understanding of the nucleon-nucleon in- 
teraction, in particular since all well established interactions fit NN scattering data with 
approximately the same precision. Substantial progress was recently achieved by the con- 
struction of a universal low energy NN potential Vj ow k based on renormalization group 
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techniques |18[. However, like the EFT potentials Vj ow k is not covariantly formulated and 
has therefore not been used in relativistic nuclear structure calculations. 

The present work applies projection techniques to map the various potentials on the op- 
erator basis of relativistic field theory which is given by the Clifford algebra in Dirac space. 
This allows to identify the different Lorentz components of the interaction and to calculate 
the relativistic self-energy operator in matter. The basic results of our approach have been 
presented in A more extended investigation is in preparation. The philosophy behind 
this approach is based on the fact, that any NN interaction, whether relativistic or not, 
is based on a common spin-isospin operator structure which invokes certain scales: long- 
range spin-isospin dependent forces, essentially given by the one-pion exchange, short- and 
intermediate-range spin-independent interactions, short-range isoscalar spin-orbit interac- 
tions and quadratic isovector spin-orbit interactions. 



II. CO VARIANT AMPLITUDES 

In a covariant formulation the spin-isospin structure of the interaction is constrained by 
the symmetries of the Lorentz group. E.g. the Born scattering matrix of covariant OBE 
potentials is given by the sum over the corresponding scalar, pseudoscalar and vector mesons 
a 

V(q',q)= ^ D a tf ~ q) T$ , (1) 

a=s,ps,v 

where g M and q'^ are the cm. momenta of the incoming and outgoing nucleons. T a are form 
factors applied to the vertices T a . The meson propagators read 

°~tf-«) = V-«)'-<. P) 

for scalar and pseudoscalar mesons s,ps and vector mesons v. Isospin factors (ri-r 2 ) are 
suppressed in The Dirac structure of the potential is contained in the meson- nucleon 
vertices 



T s = g s l, T ps = gps^f'ilb, r„ = g v Y + . (4) 

For the pseudoscalar mesons tt and rj a pseudovector coupling is used and the u meson has 
no tensor coupling, i.e. /« = 0. The potential V(q', q), i.e. the OBE Feynman amplitudes 
are obtained by sandwiching V between the incoming and outgoing Dirac spinors. 
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From a low energy expansion of the full field-theoretical OBE Feynman amplitudes into a 
set of spin and isospin operators one obtains the representation of non-covariant potentials 

v=Yy i +v{T V T 7 \o i . (5) 

i 

where the operators Oi, assuming identicle particle scattering and charge independence, are 
given by 

01 = l, 

02 = cr l' cr 2? 

3 = (ovk)(<7 2 .k), (6) 

4 = !(<ri + <t 2 ) • n, 

5 = (cr r n)(<x 2 -n), 

where k = q' — q, n = qxq' = Pxk and P = |(q + q') is the average momentum. 
The potential forms are then functions of k, P, n and the energy. In order to perform a 
non-relativistic reduction, usually the energy E is expanded in k 2 and P 2 

/W 2 \ 3 W 2 P 2 

E(q) = — + P 2 + M 2 ) ~ m + — + — . (7) 

and terms to leading order in k 2 /M 2 and P 2 /M 2 are taken into account. The meson prop- 
agators D a (q' — q) given in Equation (J2J) and (j3J) are approximated by their static form. 

The modern non-covariant Nijm 93 as well as the Nijm I and II potentials are constructed 
from approximate OBE amplitudes and thus based on the operator structure (JEJ). However, 
since they are separately fitted in each partial wave (except Nijm 93) they are often referred 
to as phenomenological models. The Argonne AVis 3 arid the Reid93 potential are purely 
phenomenological in the sense that the OBE picture is released. Based on the symmetries of 
(JHJ), the intermediate and short-range part is parameterised by phenomenological functions 
V a and only one-pion-exchange is contained explicitly. The EFT potentials, containing one- 
and multi-pion exchange explicitely, are even more rigorous since most part of the nuclear 
repulsion is carried by regularizing counter terms. Such an approach is justified by the fact 
that heavy meson (p, u) exchange cannot be resolved up to the momentum scale of about 
400 MeV which is constrained by NN scattering data. We apply the EFT Idaho potential 
which fits NN scattering data with similar quality as AVis or the Nijmegen potentials. 
The same philosophy is behind Vj ow k which can be viewed as the condensation of the various 
formulations to a model independent result. 

Although based on the same operator structure as OBEPs, it is not straightforward to fix 
the Lorentz character of the various pieces of the non-covariant interactions AVis, Reid93, 
EFT (Idaho) and V low k . 
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However, any two-body amplitude can be represented covariantly by Dirac operators 



and Lorentz invariant amplitudes 20]. A relativistic treatment invokes automatically the 
excitation of anti-nucleons. However, NN scattering, in both, relativistic and non-relativistic 
approaches is restricted to the positive energy sector and neglects the coupling to anti- 
nucleons. As a consequence one has to work in a subspace of the full Dirac space where 
on-shell two-body matrix elements can be expanded into five Lorentz invariants. A possible 
choice of a set of five linearly independent covariant operators are the scalar, vector, tensor, 
axial-vector and pseudo-scalar Fermi covariants 

r m = {S,V,T,P,A} (8) 

with 

S = 1(8)1, V = 7^® 7At , T = <7^<8><v, P = 7s®75, A = 7 57 m ®757m- (9) 

Working with physical, i.e. antisymmetrized matrix elements V A , one has to keep in mind 
that direct and exchange covariants T m (where Dirac indices of particles 1 and 2 are inter- 
changed) are coupled by a Fierz transformation. Since the low momentum part of the NN 
interaction is totally dominated by the pseudovector one-pion exchange (OPE) it is prefer- 
able to choose an operator basis where the pseudovector exchange is completely separated 
from the remaining operator structure and allows thus a more transparent investigation of 
the short- and intermediate-range parts of the potentials which are actually the interesting 
ones. 

This can be achieved by the following set of covariants originally proposed by Tjon and 
Wallace 

r m = {S, -S, (A - A), PV, -PV} . (10) 

PV and PV are the direct and exchange pseudovector covariants, analogous to the pseu- 
doscalar covariant P, however, with 75 replaced by (jf — ^)/2M 7 5. Thus the on-shell 
(|q| = |q'|) scattering matrix is given by 

V I {c(^) = Y,9 I m {H,0)T m , (11) 

m 

where 9 is the cm. scattering angle and I = 0, 1 the isospin channel. For the Hartree- 
Fock self-energy it is sufficient to consider 6 = when antisymmetrized matrix elements 
are used since 9 = it contains then only redundant information. The transformation of the 
Born amplitudes from an angular-momentum basis of a given NN potential, where the 1 Sq 
partial wave of the studied potentials is shown in Fig. ^ onto the covariant basis (jlljl is now 
standard and runs over the following steps: 

\LSJ) — > partial wave helicity states — > plane wave helicity states — > covariant basis. 
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FIG. 2: Isospin-averaged Lorentz invariant amplitudes cj^(|p|,# = 0) for the different NN po- 
tentials after projection on the Dirac operator structure. The pseudovector representation of the 
relativistic operator basis is used. As a reference the amplitudes from solely OPE and from a + u 
exchange, both with Bonn A parameters, are shown. 

The first two transformation can e.g. be found in Refs. jEjj. The last step has to be 
performed numerically by matrix inversion 0, E3 ■ 

Fig. |21 shows the resulting isospin-averaged amplitudes g m = ^(g^° + •ig^ 1 ) at 9 = for 
the various potentials. Since these amplitudes are not very transparent quantities, Fig. 121 
includes as a reference in addition the contributions from only OPE and from only a and u 
exchange, both taken from Bonn A. The four amplitudes g$, g$, g^ and g^y are very close 
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for the OBEPs Bonn A, Nijm 93, CD-Bonn, Nijmegen I and II, and the non-relativistic 
AVis, Reid93, Idaho and Vj ow k- The direct pseudovector amplitude gpy falls somewhat out 
of systematics. This amplitude is, however, of minor importance since it does not contribute 
to the Hartree-Fock self-energy (fT^j) and to the single particle potential. The dominance 
of the OPE at low |q| is reflected in the pseudovector exchange amplitude g^y which is at 
small |q| almost two orders of magnitude larger than the other amplitudes. 

For the covariant OBEPs Bonn and CD-Bonn and the OBE-like Nijmegen potentials a 
coincidence at the level of covariant amplitudes is not completely unexpected. However, it 
is remarkable that this agreement transfers to the EFT Idaho potential and also to Vj ow k 
which are of completely different character and theoretical background. 

The high momentum part of the interaction, on the other hand, is dominated by heavy 
meson exchange in the OBEPs and the corresponding amplitudes g$, ^S' 9 A approach the 
a + uj exchange result. This is also true for AVi 8 and Reid93 which could not have been 
expected a priori. Deviations from the a + uj amplitudes, e.g. due to exchange of isovector 
mesons p and 5 in the OBEPs or the corresponding isovector operators in AVig are moderate 
at large |q|. These deviations are more pronounced at small |q|. In summary we find a 
remarkable agreement between the OBE amplitudes and those derived from AVis or Reid93. 
This means that for on-shell scattering AVig can be mapped on the relativistic operator 
structure where the local phenomenological functions Vi, Eq. (jSJ), play the same role as the 
meson propagators plus corresponding form factors in the meson exchange picture. 

Turning now to the chiral Idaho potential, where the momentum-space NN amplitude is 
also based on the operator structure given in Eq. El the situation is different. At low |q| 
the amplitudes derived from the Idaho potential behave qualitatively and quantitatively like 
the previous ones, i.e., they are very close to Bonn A, CD-Bonn and AVis- Again we would 
conclude that also the chiral potential can be mapped on a relativistic operator structure. 
The functions Vi and V[ in combination with the corresponding operators, derived from 
fourth order 2n exchange plus contact terms, lead to a structure which is similar to that one 
imposed by the OBE picture. However, clear deviations appear in the cut-off region between 
400 and 500 MeV. The short-range interactions are strongly suppressed by the exponential 
cut-off form factors and as a consequence the Idaho approaches rapidly the OPE result for 
momenta above 400 MeV. The Vj ow k potential is only displayed up to a momentum of 400 
MeV since the high momentum nodes are integrated out down to a scale of A = 2.1/m~ 1 
in order to arrive at a low-momentum potential. It shows a similar behaviour compared 
to the other potential models especially when compared to the chiral EFT potential. As 
in the chiral EFT potential where the short-range physics is absorbed in contact terms, 
the unresolved short-range part is removed and replaced by contact interactions. Therefore 
when going to higher |q| the amplitudes of the Vj ow k potential deviate from the a + uj 
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exchange. Nevertheless the V\ ow k matrix elements are provided in numerical rather than 
analytical form. When mapped on a general relativistic operator basis one finds that the 
Dirac structure is comparable to that of the chiral EFT potential, the relativistic and non- 
relativistic OBEPs and phenomenological potentials. 



III. SELF-ENERGY 



With the covariant amplitudes at hand, one is now able to determine the relativistic 
mean field in nuclear matter by calculating the relativistic self-energy E in Hartree-Fock 
approximation at tree level. We are thereby not aiming for a realistic description of nuclear 
matter saturation properties which would require a self-consistent scheme. Moreover, short- 
range correlations require to base such calculations on the in-medium T-matrix rather than 
the bare potential V |23j . This leads to the relativistic Dirac-Brueckner-Hartree-Fock scheme 
which has been proven to describe nuclear saturation with quantitatively satisfying accuracy 

The self-energy for the nucleon k follows from the interaction matrix V by integrating 
over the occupied states q in the Fermi sea 



E a/3 (/c, k F ) 



(27T) 



(12) 



The Dirac propagator 

G D (q) = [4 + M]2m5(q 2 - M 2 )Q{q )Q(k F - |q|) (13) 

describes the on-shell propagation of a nucleon with momentum q and energy = 
^/q 2 + M 2 inside the Fermi sea. In isospin saturated nuclear matter the self-energy consists 
of a scalar E s , a time-like vector E Q and a spatial vector part E v 



E(fc, k F ) = E s (/c, k F ) - 7o E Q (fc, k F ) + 7 • kE v (fc, k F ), 



(14) 



which are given by 
E a 



kF d 3 q M 



4#s -9§ + 4 5 'a - 



(k* 1 - q») 



M \2 



(2tt) 3 
kF d 3 q k-q 
(2vr) 3 |k| 2 £ q 



E k (k" - q 
AM 2 



AM 2 

m 



-9vv 



(15) 



k z {MP - q") 2 



Fig. El (a) shows the tree level self-energy components obtained with the various potentials 
in nuclear matter at saturation density with a Fermi momentum k F = 1.35 fm _1 . All 
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potentials yield scalar and vector fields S s and S Q of comparable strength: a large and 
attractive scalar field S s ~ —(450 -r- 400) MeV and a repulsive vector field of S Q ~ +(350 4- 
400) MeV. These values are comparable to those derived from RMF phenomenologically 
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FIG. 3: Tree level results for (a) the scalar S v and vector S Q self-energy components and (b) 
spatial vector self-energy component kS v in nuclear matter at kp = 1-35 fm for the different NN 
interaction models. 



and also from QCD sum rules. Also the explicit momentum dependence of the self-energy 
is similar for the various potentials. The Idaho mean fields follow the other approaches at 
low k but show a stronger decrease above k ~ 2 fm -1 which reflects again the influence of 
the cut-off parameter. The spatial components kS v are shown in Fig. |3] (b). Also here the 
various potentials agree quite well. 

As known from self-consistent DBHF calculations the spatial vector self-energy is 

a moderate correction to the large scalar and time-like vector components S s and S Q . This 
is already the case at tree level where kS v is about one order of magnitude smaller than 
the other two components. The spatial self-energy originates exclusively from exchange 
contributions, i.e., the Fock term, and vanishes ,e.g., in the mean field approximation of 
RMF theory. 
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FIG. 4: Single particle potential in nuclear matter at k-p = 1.35 fm _1 , determined from the tree 
level Born amplitudes of the various potentials. The single particle potential determined from the 
relativistic self-energy components after projection onto the covariant operator basis is compared 
to a non-relativistic calculation (stars) where partial wave amplitudes are summed up directly. 

In Fig. E| the single particle potential is shown, which is defined as the expectation value 
of the self-energy E 

< u0)|7°£|«0) > M 



U B . p .(k) 



<u{k)\u{k)> E(k) 



< u(k)\E\u(k) > 



and reads 



M 

U s . p .(k,k F ) = — S s 



k^ 



MS, 



S v k 2 



(16) 



(17) 



E Vk 2 + M 2 Vk 2 + M 2 

It reflects the well known fact that various two-body potentials are rather different al- 
though they are phase-shift equivalent, i.e. they describe NN scattering data with about the 
same accuracy when iterated in the Lippmann-Schwinger equation [23]. The differences at 
tree level are mainly due to differences in the short-range part of the interaction [IS]. Fig. 0] 
includes also the results from a 'non-relativistic' calculation of U s . p . where the partial wave 
amplitudes are directly summed up. In particular at low momenta the non-relativistic and 
the relativistic calculations show an excellent agreement which demonstrates the accuracy 
of the applied projection techniques. 



11 



In order to estimate the influence of short-range correlations and self-consistency, Fig. 
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FIG. 5: Tree level scalar and vector self-energy components in nuclear matter at = 1-35 fm -1 
are compared to corresponding values from a full self-consistent relativistic Brueckner (DBHF) 
calculation. In both cases the Bonn A potential is used. 

|S] compares the tree level result from Fig. |3] for Bonn A to a corresponding full DBHF 
calculation at hp = 1.35 fin -1 . For DBHF the approach of is used (subtracted T-matrix 
in pv representation). The DBHF calculation yields reasonable saturation properties with a 
binding energy of -Ebind — —15.72 MeV and a saturation density of p = 1.84 fm~ 3 [7]. It is 
no doubt that higher order correlations are essential for saturation of nuclear matter. The 
correlations lead to a general reduction of the vector self-energy by a shift of about 50 MeV. 
Self-consistency and correlations also weakens the momentum dependence, in particular for 
S s . However, except of the 50 MeV shift of S , the absolute magnitude of the self-energies 
is not strongly modified in the realistic calculation. This means that one can expect that 
the large attractive scalar and repulsive vector mean fields will also persist for the other 
interactions when short-range correlations are accounted for in a full relativistic many-body 
calculation. 




12 



IV. SUMMARY 



We presented a model independent study of the Dirac structure of the nucleon-nucleon 
interaction. The potentials were projected on a relativistic operator basis in Dirac space 
using standard projection techniques which transform from a partial wave basis, where the 
potentials are originally given, to the basis of covariant amplitudes. The different approaches 
can now be compared at the level of these covariant amplitudes and, moreover, this allows 
to calculate the relativistic self-energy operator in nuclear matter. For both, the covariant 
amplitudes and the tree-level Hartree-Fock self-energy, we observe a remarkable agreement 
between the OBEPs (Bonn, CD-Bonn, Nijmegen), the phenomenological AVis and Reid93 
potentials and the EFT based Idaho and V\ ow k potentials. The structure of the interaction 
enforces the existence of large scalar and vector fields as a model independent fact. The scale 
of these fields is set at the tree level. Although essential for nuclear binding and saturation, 
higher order correlations, i.e. Brueckner ladder correlations, change the size of the fields by 
less than 20% 0]. The magnitude of the fields is similar to that predicted by relativistic 
mean field phenomenology, relativistic many body correlations and also by QCD sum rules. 
We conclude that the appearance of large scalar and vector fields in matter is a general and 
model independent consequence of the vacuum structure of the NN interaction. Relativistic 
dynamics is therefore essential for nuclear systems. 
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